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Resume. Dans le contexte de la propagation d'ondes elastiques dans des milieux endom- 
mages, une approche analytique pour la diffusion d'ondes antiplanes par des reseaux periodiques 
bidimensionnels de fissures est developpee. Avant d'envisager I'etude de reseaux de fissures, la dif- 
fusion simple d'une onde elastique antiplane par une fissure plate est prealablement etudiee. En- 
suite, en employant la formule de representation d'une fissure plate et la condition de periodicite 
d'une rangee de fissures situees sur la meme droite antiplane, une equation integrale de frontieres 
est alors obtenue pour le pour le saut de deplacement a travers une fissure dite de reference. 
Des resultats numeriques pour les coefficients de refiexion et de transmission sont presentes en 
fonction de I'espacement entre fissures, de la frequence d'excitation et de Tangle d'incidence. Fi- 
nalement, I'etude de la propagation d'une onde antiplane au travers d'un empilement periodique 
d'un nombre fini de reseaux lineaires (i.e. un reseau dit bidimensionnel periodique) est axee sur 
I'analyse des bandes spectrales interdites et permises. Les effets dus a une perturbation de la 
periodicite ont egalement ete analyses. 

Abstract. In the context of elastic wave propagation in damaged solids, an analytical ap- 
proach for scattering of antiplane waves by two-dimensional periodic arrays of cracks is developed. 
Before considering the study of arrays of cracks, the scattering of an antiplane wave by a flat 
crack is first studied. Then, using the representation formula for the scattered displacement by 
a flat and by considering the periodicity condition of the crack-spacing, a boundary integral 
equation is obtained for the crack face displacement of the reference crack. Numerical results for 
the reflection and transmission coefficients are presented as functions of the crack-spacing, the 
frequency of excitation, and the angle of incidence. Finally, the propagation of antiplane waves 
by two-dimensional periodic arrays of cracks is studied. Despite the use of a finite number of 
linear arrays, one recognizes the effects of band-pass filtering or band rejection characteristics of 
the transmission spectra of a periodic medium. Effects due to a disorder in the periodicity are 
also analysed. 

Mots-cles : endommagement ; reseau de fissures ; propagation d'ondes elastiques 
Keywords : damage ; array of cracks ; elastic wave propagation 
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Abridged English version 



Localized damage in form of distributed coplanar micro-cracks often occurs at the bounding 
surfaces between two materials of the same or different mechanical properties. Such micro-cracks 
may be induced by manufacturing, materials processing, or under service conditions. The nu- 
cleation, growth and coalescence of that localized damage may gave rise to substantial stiffness 
and strength reductions of the materials. Thus, the detection and characterization of the damage 
state in an engineering structure is of considerable interests for improving the manufacturing 
and the material processing, estimating the life-time, and avoiding catastrophic failures of a 
damaged structure. This research has also considerable interest, on the one hand to ultrasonic 
quantitative non-destructive evaluation for detecting and characterizing the damaged state of 
composites, and on the other hand to linear elastic fracture mechanics to assess the initiation 
and growth of a developed macro-crack under dynamic loading conditions. With this motiva- 
tion, the present analysis is concerned with elastic antiplane (shear horizontal) wave propagation 
in isotropic materials containing micro-cracks arranged periodically in two space-directions. Al- 
though an antiplane wave has limited applications in ultrasonic non-destructive evaluation, some 
useful information can be extracted from this simple model. 

Summary 

A method similar to that proposed by Zhang and Gross [Tl is used to formulate the problem 
of scattering by a periodic array of collinear shear cracks. The geometry of the problem is 
presented in Fig. 1. By using a representation formula for the scattered displacement and by 
considering the periodicity condition of the crack-spacing, an integral equation is obtained for 
the crack face displacement of a reference crack, which is solved numerically by using a Galerkin 
method. In general, the scattered wave field consists of an infinite number of wave modes, each 
with its own cut-off frequency. Under the cut-off frequency the wave modes are standing wave 
modes which decay exponentially with distance to the array. At frequencies higher than the cut- 
off frequency the wave modes will be propagating. At the cut-off frequency a mode conversion 
from a standing into a propagating wave mode occurs. Sufficiently far from the array, and at 
sufficiently low frequencies, only the zeroth (lowest) order mode is dominant, and that mode 
represents homogeneous plane waves corresponding to reflected and transmitted waves by the 
cracks. Numerical results for the reflection (transmission) coefficient are presented, as function 
of the crack-spacing, the wave frequency, and the angle of incidence. 

The wave propagation in two-dimensional periodic arrays of flat cracks is studied next. The 
global array of Fig. 5 is composed of an arbitrary number S of infinite periodic arrays of collinear 
cracks, regularly spaced by a distance D along the axis 2/3. Its thickness is h = {S — \)D. We 
assume that all linear arrays are identical, with d their common period (along the axis yi) and 2a 
the crack width. The reflection and transmission coefficients of the linear array s (s = 1, 2, S) 
are known from the previous investigation. From these coefficients, the global reflection and 
transmission coefficients of the two-dimensional array are calculated without difficulty by using 
an iterative algorithm. Other closely related problems, some more complicated as they involve 
inplane ((longitudinal or shear vertical)) waves, can also be found in the literature; see, e.g. 
[2] and [3]. The formalism outlined in the text is valid only at low frequencies and only the 
lowest order modes can propagate in the array. The other modes are evanescent and decrease 
exponentially with the propagation distance. They can therefore be neglected if the distance D 
is large enough, such that an evanescent mode can not disturb significantly a neighbor linear 
array. Despite the use of a finite number of linear arrays, one recognizes the effects of band-pass 
filtering or band rejection characteristics of the transmission spectra of a periodic medium. The 
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band rejections occur in the vicinity of frequencies for which there is coincidence between the 
spacing of hnear arrays and muUiples of the half-wavelength. 

Finally, for an infinite periodic two-dimensional array of cracks, the effective wave number of 
antiplane waves propagating in the homogeneous medium equivalent to two-dimensional array is 
calculated by using the Floquet theory. 

The theory developed here can be further generalized to a two-dimensional array of cracks 
with any type of periodicity, such as a variation of the crack-sizes or a random distribution of 
linear arrays. 

1 Introduction 

Un endommagcmcnt localise sous forme d'une distribution de microfissures coplanaires se 
produit souvent au niveau des surfaces de liaison entre deux materiaux aux proprietes mecaniques 
identiques ou differentes. De telles microfissures peuvent etre induites durant la fabrication, par 
les materiaux traitants, ou bien les conditions d'utilisation. La germination, la croissance et la 
coalescence des dommages localises peuvent provoquer des reductions substantielles de rigidite et 
de resistance des materiaux. Ainsi, la detection et la caracterisation de I'etat d'endommagement 
dans une structure d'ingenierie est particulierement importante pour ameliorer le processus de 
fabrication, les materiaux traitants, et eviter les ruptures catastrophiques. Cette recherche a 
aussi un interet considerable, d'une part pour revaluation ultrasonore non-destructive, dans le 
but de detecter et caracteriser I'etat d'endommagement de composites, et d'autre part, pour la 
mecanique elastique lineaire de la rupture qui permet d'evaluer le declenchement et la croissance 
d'une « macro-fissure > pouvant se developper dans des conditions de charges dynamiques. Bien 
que beaucoup de facteurs puissent surgir in situ, les composites unidirectionnels renforces ainsi 
que les laminees a plis croises peuvent etre commodement approximes comme des materiaux 
elastiques a symetrie isotrope ou orthotrope. Une telle approximation est possible lorsque la 
longueur d'onde de I'onde incidente est beaucoup plus grande que les longueurs caracteristiques 
des materiaux composites (espacement, diametre des fibres dans un composite renforce, epaisseur 
de differentes couches d'un stratifie). 

Ainsi, I'analyse de la section 2 concerne la propagation d'une onde antiplane dans des materiaux 
isotropes contenant des microfissures disposees periodiquement dans une direction spatiale, et 
formant un reseau lineaire infini. La section 3 a pour objectif I'etude de la propagation d'ondes 
antiplanes dans des reseaux periodiques a deux dimensions spatiales. La methode de calcul des 
coefficients de reflexion et de transmission du reseau bidimensionnel, ainsi que I'obtention des 
courbes de dispersion et d'attenuation correspondantes, reposent sur la decomposition des reseaux 
bidimensionnels en un nombre flni de reseaux lineaires dont leurs propres coefficients de reflexion 
et de transmission sont desormais connus. 

2 Arrangement periodique de fissures colineaires 

Dans cette section, un modele est presente afin de calculer les coefficients de reflexion et 
de transmission par un arrangement periodique de fissures colineaires dans un solide elastique 
isotrope, pour le cas d'une onde antiplane. Meme si les ondes antiplanes ont des applications 
limitees dans revaluation ultrasonore non-destructive, des informations utiles peuvent etre ex- 
traites a partir de ce modele simple. Une methode similaire a celle proposee par Zhang et Gross 
[1] est employee pour formuler le probleme de diffusion par un reseau de fissures lineaire et in- 
fini. Les deplacements diffuses sont exprimes sous forme d'integrales de Fourier contenant les 
deplacements au niveau des levres des fissures. En employant la formule de representation d'une 
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fissure plate et la condition de periodicite d'une rangee de fissures situees sur la meme droite 
antiplane, une equation integrale de frontieres est alors obtenue pour le deplacement des levres de 
la fissure dc reference. Gcncralement, le champ de deplacement diffuse se compose d'un nombre 
infini de modes d'onde, chacun avec sa propre frequence de coupure. Au-dessous de leur frequence 
de coupure, les modes sont dits ondes stationnaires et decroissent exponentiellement au cours de 
la propagation. Au-dcla de la frequence de coupure, les modes deviennent des ondes homogcncs. A 
la frequence de coupure, une conversion du mode stationnaire en un mode homogene se produit. 
Lorsque situe sufHsamment loin de I'interface, et a des frequences sufRsamment basses, le mode 
d'ordrc zero est dominant, reprcsentant les ondes planes homogcncs rcflcchies et transmises par 
les fissures. Des resultats numeriques pour les coefficients de reflexion et de transmission sont 
presentes en fonction de I'espacement entre fissures, de la frequence d'excitation et de I'angle 
d'incidence. 



2.1 Champ diffuse par une fissure 

Considerons un solide elastique isotrope, homogene, infini, contenant une fissure de longueur 
2a. La fissure est centree a I'origine du systeme (j/i, et est infiniment longue dans la direction 
perpendiculaire y2 ■ Nous soumettons cctte fissure a un champ acoustique antiplan, defini par un 
champ de deplacement dans la direction y2., de la forme : 

uT (yuys) = uoe"=^(2'i'''"^»+«^'=°''^»), (1) 

ou uo et 00 representent respectivement I'amplitude de I'onde et Tangle d'incidence. D'apres la 
loi de Hooke, la contrainte 0-23'^ associee est donnee, dans le plan de la fissure {ys = 0), par 

a^f (yi , 0) = iuokTfXo cos Ooc'^-y^ '^^ , (2) 

oil kr = lu/ct est le nombre d'onde, ct est la celerite transverse et /xq est le module de rigidite 
du solide. Puisque il y a une invariance en la translation dans la direction 2/2, le probleme de 
diffusion resultant dcvient natiirellement bidimensionnel. 

Quand la fissure interfere avec I'onde incidente, des ondes diffusees sont generees. Ces ondes 
produisent un deplacement diffuse antiplan •U2'*) solution de I'equation de Helmholtz scalaire 



iVuVs) + (2/1,2/3) + kluf {yuya) = 0. (3) 

Notons que le deplacement diffuse est antisymetrique par rapport au plan (t/i, ^2)- Nous pouvons 
done limiter I'etude au demi-espace ys > 0. En supposant que les levres de la fissure sont fibres 
de contraintes, la condition aux limites du probleme se traduit par : 

'^21' (2/1,0+) ^ -iuofcT/iocos0oc"^-^^^"^«, (I2/1I < a). (4) 
En dehors de la fissure, le champ de deplacement vaut : 



.dif 



(2/1,0+) = 0, (|2/i|>a). (5) 



2.1.1 Representation integrale 



de uf^ : 



Definissons la transformee de Fourier spatiale (par rapport a la variable yi) directe et inverse 

uf (e, 2/3) = ^ / «f (y) e'«^Myi, (6a) 
v27r J-00 
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[y] - -= / uf (e, ys) e-'«^Me (6b) 
Apres avoir applique la transformee (|6a| a I'Eq. ([3|) et defini la fonction /3 telle que : 



^{ki-eY'\ e<ki, 

Tequation de Helmholtz ([3]) devient une equation differentielle ordinaire de variable j/3 



^(e,2/3)-/?'(e)"f (^,2/3)=0. (8) 
Une solution bornee dans la region j/3 > et satisfaisant la condition de Sommerfeld [5], s'ecrit : 

uf{ty3)^CiOe-^y^, (9) 
oil la fonction C est donnee par I'Eq. (pa|) sous la forme : 

C (0 = r uf (yi, 0+) e'«^Myi. (10) 

Nous pouvons maintenant substituer I'Eq. (flUl) dans I'Eq. © et appliquer la transformee inverse 
(j6b|) . La solution du probleme de diffusion (H))-® pent s'ecrire sous la forme : 

uf (yi, 2/3) = ^sgn (yg) / / uf {v, 0+) c'«''-^^)e-^(«)l^^ldz;de (H) 

L'equation (jlip montre qu'il est possible d'evaluer wf'^ quel que soit (yi, ya) e M^, si Uj'^ (w, 0+) 
est connu pour \v\ < a. En utilisant I'Eq. pip dans la condition aux limites nous obtenons 
l'equation integrale suivante pour Wj'* (w, 0^) : 



(«,0+) / /3(e)e'«(^~2^iMCdw = fcTCOse'oe''^-^yi"'"^«, (|yi| < a). (12) 



Angel [6j a montre que I'integrale double dans le membre de gauche de I'Eq. (fT2)) se reduit 
a une integrale singuliere de type Cauchy. L'equation integrale ((T2|) ne pent pas etre resolue 
analytiquement, mais il est possible de la resoudre numeriquement en utilisant, par exemple, une 
formule d'integration de type Gauss-Tchebychev (voir '7], chapitre 2, pour plus de details). 



2.1.2 Champ lointain 

Soit (r, 6) les coordonnees polaires deflnies par — yl+yl et cos 9 = ys/r. En champ lointain, 
lorsque r 3> a, le deplacement diffuse (1111) pent s'ecrire comnie le produit d'une fonction radiale, 
qui depend de kxr, et d'une fonction de I'angle d'observation 0, qui depend a priori de Tangle 
d'incidence Oq et de la frequence. Nous avons [TD] : 

uf (r, 9, 9o) = uo^f^e^'^^f {9, 9^) , (r » a, 9e [0, 2^)) , (13) 
V KTr 

oil V amplitude de diffusion en champ lointain f est definie par : 

f{9,9o)^kTCot9 f (w,0+)e-"=^'"'''^^dv. (14) 
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2.2 Champ diffuse par un arrangement periodique de fissures co- 
fineaires 

L'objectif de cette section est d'ecrire le champ diffuse par un reseau lineaire infini sous la 
forme d'ondes planes reflechies et transmises. La geometric du reseau et les notations utilisees 
sont presentees sur la Fig. 1. 

2.2.1 Enonce du probleme 

Le reseau est compose d'une disposition periodique de fissures antiplanes colineaires. Soit 2a 
la longueur de chaque fissure et d la distance entre les centres de deux fissures voisines (pas de 
reseau). Soit (CiiCs) le repere « local » lie a une fissure. Le repere « global » (yijj/s) est relie 
au repere local (Ci, Ca) P^r : 

D'apres I'Eq. (fTTj) . le champ de deplacement difi^use au point d'observation (2/1,2/3) est done 
donne par : 

uf{yi,y3)^^sgn{y,) J2 / / ^^UCi,0+)e'«(^^+^'^"^^)e-^(«)l«='ldCide, (16) 

Oil U2 (Ci, 0+) est le deplacement diffuse sur la fissure j. 

2.2.2 Hypothese de periodicite 

D'une fissure a une autre le deplacement incident est le meme mais dephase de la quantite 
/cT^sin^o- Ce terme vaut 1 quand 9o — 0, c'est-a-dire en incidence normale, quand toutes les 
fissures sont excitees en meme temps. En substituant (TTSl) dans (0), le champ de contrainte 
incident sur les levres des fissures pent s'ecrire comme : 

ivi , 0) = ^f^ (Ci ) e"^-^'^ ^" , (17) 
oil (T23 (Ci) est le champ de contrainte incident sur la fissure de reference {\Ci \ < a, Ca — 0) 

'^23 (Ci) = iuokrfio cos^oe''^-^^ (18) 

Par suite, pour respecter I'hypothese de periodicite dans la direction j/i, nous pouvons ecrire 
le deplacement U3 (^1,0+) sur la forme [8^ : 

^.^2(Ci,0+) = U2(Ci)e"=-^'*^'"«°, (19) 

oil U2 (Ci) (|Ci| < o) est le deplacement diffuse sur la fissure de reference et j = 0,±1, .... De 
cette fagon, nous pouvons eliminer un grand nombre d'inconnues de I'Eq. en ne conservant 
que le deplacement diffuse sur la fissure de reference. 
En substituant ([T^ dans nous obtenons : 

«f (yi,2/3) = ^sgn(y3) / / ^2(Ci)e'(*^-^"^''+«)^'^e'€(^^-^^)e-^(«)l^3ld(;^de (20) 
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Rappelons maintenant les proprietes de la fonction de Dirac 6 : 

oo oo 

277 <5[(fcTsin0 + Od-2j7r] = ^ ^iKkT sin 4>+i)d . ^21) 

J — — OO j— — oo 

S{x) = \d\S{dx), Vx, yd^ const eC] (22) 

^(C-eo)<P (Ode = V^GCi(C). (23) 

Nous obtenons alors, compte tenu de I'Eq. ((2T|) : 

oo /.oo /.a 



'f(2/i,y3)=sgn(y3) ^ / / W2 (Ci) -5 [(fcr sin^o + - 2i^] e'«(^^"^^)e-'^(«)l^^'ldCide 

(24) 



Soit, en utilisant I'Eq. ([22]) : 



1 ^ /"OO 

f (2/i,y3) = ^sgn(y3) ^ y_ <5(e-«,)^(e)de, (25) 



j = -oo 



avec 



(p{0= U2(Ci)e'«('^i-^i)e-''(«)l2'-'ldCi (26) 

J —a 

et 

ajd = 2j7r — fcrdsin^o- (27) 

Enfin, il resulte des Eqs. (P5|) - (P7)) apres avoir utilise I'Eq. que le champ de deplacement 

diffuse par I'ensenible des fissures est de la forme : 

4''(yi,2/3) = Jsgn(2/3)y U2{Ci)S{Ci-yi,y3)dCu (28) 

avec 

oo 

S{x,y)= J2 e'"^-^-^^l^l, (29) 

j=-oo 

oil 

Pj=P{a,). (30) 

L'equation ((28)) montre qu'il est possible d'evaluer Uj'^ en tous points de I'espace, si le champ 
diffuse U2 (Ci) est connu pour |Ci| < a. Par consequent, U2 sert de quantite inconnue fondamen- 
tale. 

Le champ de contrainte associe au champ de deplacement (|28p s'ccrit : 

^2^ (yi,2/3) = ^ r U2 (Cl) S' (Cl - yi,2/3)dCl, (31) 
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avec 



S'{x,y)^- p^e^'.-PM, (32) 

En faisant s'approcher le point d'observation {yi, y^) des levres de la fissure de reference dans 
I'Eq. (|31|) . et en utilisant la condition aux limites (|4|), nous obtenons Tequation integrale suivante 
pour Finconnue U2 : 

/a oo 
"2(Ci) 5^ /3,e'"^^^^""^^dCi =-iuofcTdcos0oe"=-^^^'"'°, \yi\ < a. (33) 

i=-oo 

L 'equation integrale (|33p pent etre resolue numeriquement en adoptant une methode dite de 
Galerkin (voir [T] pour plus de details). La solution M2 de I'equation integrale est ainsi 
determinee par un calcul exact de diffusion multiple prenant en compte toutes les interactions 
possibles entre les fissures. 

2.2.3 CoefRcients de reflexion et de transmission du mode fondamental 



Dans I'Eq. (I28p . chaque j correspond a un mode d'onde. Ceci implique que le champ de 
deplacement diffuse consiste en un nombre infini de modes, chacun avec sa propre frequence de 
coupure. Pour a| > oii : 



— ; kt sm 0^) 

a 



> kT, (34) 



les modes diffuses de wf'^ sont des ondes stationnaires (ou evanescentes) qui decroissent expo- 
nentiellement avec I'augmentation de la distance du reseau, tandis que pour Oj > k^, ces modes 
representent des ondes homogenes se propageant dans les directions j/3 positives (ou negatives). 
Aux frequences de coupure aj = k^, un mode stationnaire se converti dans un mode homogene. 
On pent montrer que, suffisamment loin du reseau et a des frequences suffisamment basses, le 
mode d'ordre zero (j = 0) sera dominant. Ces ondes d'ordre zero sont les ondes planes ho- 
mogenes reflechies et transmises, qui se propagent dans la meme direction que I'onde incidente. 
Ceci devient clair en reecrivant le mode d'ordre zero sous la forme : 



(2/1, ys) = Roe'''^(v^-"^<^o~yscoseo}^ < (35) 



"f (yi,2/3) = (To-l)e''=-(2'i-"^«+^-'™^^«), y3>0, (36) 

oia i?o et Tq designent les coefficients de reflexion et de transmission du mode fondamental (d'ordre 
zero) qui sont donnes par : 

i?o = f U2 {v)e-'^^-'"''<'dv, To = 1 - i?o, (37) 

Notons que I'Eq. (1571) satisfait le bilan d'energie bien connu [9 : 

|i?o|' + iTol' = 1. (38) 
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2.3 Resultats numeriques et discussion 



Le coefRcient de reflexion du mode fondamental est trace ici sur les Figs. 2-4 comme etant 
une fonction de respacement entre fissures d = d/a, de la frequence d'excitation u} — kxa et de 
Tangle d'incidence 9o. 

La Fig. 2 montre le module du coefficient de reflexion Rq en fonction de la frequence adi- 
mensionnelle w pour trois valeurs de d et pour deux angles d'incidence (6*0 = et 6*0 = 7r/4). Le 
coefficient de reflexion est zero a w = (ondes infiniment longues), et augmente monotonement 
avec I'augmentation de a), jusqu'a ce qu'un maximum soit atteint et pour lequel les courbes af- 
fichent une pente discontinue. Cette discontinuite ainsi que les suivantes qui apparaissent lorsque 
Cj augmente, se produisent aux frequences de coupure Cj^^\ ce qui resulte des Eqs. ([7]), ^7} et 

c^a) ^ ^ _,^a)sin0o, j-1,2,.... (39) 
d 

Ce sont les valeurs de uj auxquelles un mode d'ordre superieur devient un mode propagatif 
(homogene). Pour I'incidence normale {9o = 0), I'Eq. ((5^ se reduit a w^-'^ = 2jTr/d. Considerons 
par exemple la courbe pour d = 2 et 9o = 0- Dans la gamme de frequences < a) < tt, le champ 
reflechi a une distance quelconque de — est uniquement constitue du terme j ~ dans 
I'Eq. (|^ ; les termes d'ordres superieurs peuvent etre negliges, car ils sont evanescents. Dans 
les gammes de frequences tt < a) < 27r et 27r < w < in, les termes correspondant au j — 1 
et j — 2 doivent, pourtant, aussi etre inclus dans le champ total reflechi en champ lointain. 
Seul le coefficient de reflexion |i?o| du mode fondamental est represente sur la Fig. 2. Pour une 
valeur approximative d ~ 2 (fissures etroitement alignees), le coefficient de reflexion approche 
rapidement |i?o| — 1 (reflexion totale) lorsque w augmente. Les courbes de la Fig. 2b (6*0 = 7r/4) 
ont la meme apparence generale que celles de la Fig. 2a {6q = 0). Notons pourtant que, pour les 
memes valeurs de d et w, les coefficients de reflexion semblent etre plus petits pour = 7r/4 que 
pour 9o = 0, du moins, dans la gamme de frequences < cD < 7r/2. 

Pour d — 20, I'espacement entre fissures est suffisamment grand afin pour que I'on neglige 
les interactions entre les fissures, et les coefficients de reflexion montres sur la Fig. 3 approchent 
ceux d'une fissure seule : 

|/(^^-Y°l (J»l); (40) 
noj cos 9o 

f{9,9o) etant I'amplitude de diffusion en champ lointain de la fissure, Eq. (fH]) . Notons que 
f {n — 9o,9o) = —f{9o,9o), ou f {9o,9o) represente I'amplitude de diffusion vers I'avant. Ceci 
nous fournit un controle sur I'exactitude de la methode numerique presentee ici. 

Dans la Fig. 4, la variation du coefficient de reflexion en fonction de la distance adimension- 
nelle d est montree pour une incidence normale {9q = 0), et pour plusieurs valeurs de la frequence 
adimensionnelle Q. Lorsque d s'approche de 2 (fissures etroitement alignees), nous trouvons que 
|i?o| tend vers 1 (reflexion totale). Remarquons que sur la courbe correspondant a a) = 2, un pic 
se distingue lorsque d — it. En effet, au regard de I'Eq. cD — 2 represente une frequence de 
coupure pour une distance d = n. 



3 Reseau periodique bidimensionnel de fissures 

La propagation d'ondes antiplanes dans des reseaux bidimensionnels de flssures est etudiee. 
Ces reseaux sont formes d'un empilement de reseaux lineaires infinis et periodiques. Dans un 
premier temps, nous rappelons la methode de calcul des coefficients de reflexion et de transmission 



9 



de reseaux bidimensionnels finis et periodiques. Ensuite, pour des reseaux infinis et periodiques, 
le nombre d'onde des ondes antiplanes se propageant dans le milieu homogene equivalent au 
reseau bidimensionnel de fissures est calcule a I'aide de la theorie de Floquet. 

3.1 CoefRcients de reflexion et de transmission d'un reseau fini 

Considerons le reseau global de la Fig. 5. II est compose d'un nombre arbitraire S de reseaux 
lineaires infinis, regulierement espaces d'une distance D le long de I'axe ^3. Le reseau global 
est done fini et d'epaisseur h = {S — 1) D. Nous considerons tout d'abord que tons les reseaux 
plans sont identiques, d etant la periode spatiale commune (selon I'axe yi) et 2a la taille des 
fissures. Isolons ici un reseau plan s (s = 1, 2, S) du reseau global. Le probleme de diflfusion 
d'une onde antiplane par un tel reseau a ete etudie dans la section precedente. Les coefiicients 
de refiexion et de transmission du reseau lineaire infini s sont desormais connus et notes R et 
T, oil I'indice associe au mode fondamental a ete ici supprime pour alleger les ecritures. A 
partir de ces coefficients, ceux du reseaux global, TZ et 7", peuvent etre calcules en elaborant un 
schema iteratif relativement simple. La technique employee se fonde sur celle de Huang et Heckl 
[1]. Ce formalisme est valable uniquement aux basses frequences et seuls les modes d'ordre zero 
peuvent se propager dans le reseau. Les autres modes sont des modes evanescents qui decroissent 
exponentiellement au fur et a mesure de leur propagation. 

Dans la suite, nous mettons en equation les refiexions multiples et les allers-retours entre 
deux reseaux lineaires consecutifs. Le calcul des coefficients de reflexion et de transmission du 
reseau bidimensionnel s'effectue sans difficulte en utilisant un algorithme iteratif presente par la 
suite. 

3.1.1 Reflexion et transmission par deux reseaux lineaires infinis 

Isolons ici deux reseaux lineaires paralleles consecutifs r et s du reseau global, le reseau r 
etant celui soumis a I'onde incidente. 

Le champ total transmis, compte tenu des reflexions multiples entre les reseaux r et s, s'ecrit 
a I'aide des series de Debye, sous la forme suivante : 

= uoTrTs + uoTrT.RsRre''!' + uoTrTs {RsRre''f'Y + ... 

= UoTrTs (l + R.Rre''^ + {R,Rre''^f + (41) 

= uaTrT, {1 ^ RsRre'^y\ 

Le champ total transmis a ete exprime en fonction d'une serie geometrique. La convergence 
de cette serie est assuree des lors que |i?si?re"^| < 1. Dans I'Eq. (|1T|) . nous avons introduit la 
notation : 

(/) = 2fcT-Dcos6'o. (42) 

D'une maniere similaire, le champ total reflechi par les deux reseaux lineaires pent etre obtenu 
sous la forme : 

u^^uo + uoRr + uoT^Rse'''' (l - R^Rre'''')'^ . (43) 
Les coefficients de reflexion et de transmission globaux sont alors definis par : 
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3.1.2 CoefRcients de reflexion et de transmission du reseau global 



Les formules (H^ et (pS)) peuvent etre utilisees pour elaborer une methode iterative de calcul 
des coefficients de reflexion et de transmission. Nous decomposons le reseau bidimensionnel dans 
5—1 reseaux (ceux-ci forment alors un « sous-reseau » bidimensionnel caracterise par les coef- 
ficients inconnus TZs-i et Ts-i) ct dans un autre reseau lineaire (caracterise par les coefficients 
locaux RetT). Suite a cette decomposition, les coefficients TZs et 7s , du reseau global de la Fig. 
5, peuvent s'exprimer en fonction de ceux du sous-reseau S"— 1, en remplagant (i?r, Tr) et {Rs,Ts) 
dans les Eqs. (|44l) - (|45|) par {R,T) et {TZs-i,Ts-i), respectivement. Ainsi, nous construisons les 
formules iteratives suivantes : 

ns^R + / ' f,'' (46) 

Ts = I2±:l . (47) 

Notons que les coefficients de reflexion et de transmission calcules peuvent etre valides en utilisant 
directement la relation de conservation d'energie psp . II suffit de remplacer dans cette relation 
les coefficients du reseau lineaire par ceux du reseau bidimensionnel. 



3.1.3 Exemples numeriques 

Les reseaux bidimensionnels consideres dans ce paragraphe possedent une direction de periodicite 
selon I'axe yi. Pour simplifier I'analyse des resultats, les coefficients de reflexion (transmis- 
sion) sont ici calcules en incidence normale. Deux exemples numeriques sont consideres : i) une 
repartition periodique selon I'axe j/3 ; ii) une repartition aleatoire de reseaux lineaires selon I'axe 

2/3- 



Repartition periodique de reseaux lineaires. Des resultats numeriques sont obtenus dans 
le cas d'un reseau bidimensionnel compose de S* = 1, 2, 5 et 15 reseaux lineaires infinis. La Fig. 6 
niontre I'influence du nombre des reseaux lineaires S sur le coefficient de transmission en energie 
|T|^. Pour d{— d/a) = 5 et D{— D /d) = 1, la premiere frequence de coupure correspond alors a 
kxD — 2tt (voir I'Eq. (1391) '). Observons que dans une gamme de basses frequences, le nombre des 
minima et des maxima egale le nombre des rangees des reseaux lineaires moins un. La comparaison 
des courbes pour 5 = 5 et 5 = 15 illustre la formation des bandes passantes et bandes d'arret 
(ou interdites). L'existence des bandes passantes et d'arret donnent au reseau bidimensionnel les 
caracteristiques d'un flltre mecanique [111 dans lequel des ondes ayant la frequence et le nombre 
d'onde se trouvant dans les bandes d'arret sont attenuees, alors que d'autres sont transmises. 
Meme si le nombre des reseaux lineaires est flni, de caracteristiques attendues pour un milieu 
periodique infini, telles que des bandes passantes et d'arret sont reconnues. Ces aspects sont 
etudies en detail dans la section suivante. 



Repartition aleatoire de reseaux lineaires. Considerons ici un reseau bidimensionnel com- 
pose de S" = 30 reseaux lineaires infinis et disposes periodiquement avec une periode D = 1. 
Chaque reseau lineaire est caracterise par d = 5. L'objectif de ce paragraphe est d'identifier I'ef- 
fet d'un defaut de position aleatoire de chacun des reseaux lineaires par rapport a leur position 
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initiale periodique. Seuls les premiers et derniers reseaux lineaires sont supposees fixes. Nous 
introduisons une distance liniite d'approche entre les reseaux, notee 6, dans le but de quantifier 
ce defaut de periodicite. Ainsi, b = D correspond a un reseau periodique et, a I'extreme, 6 = a 
un reseau aperiodique oil les fissures de deux reseaux lineaires consecutifs sont susceptibles d'etre 
en contact. Nous nous assurons ici que b > AD pour que la contribution d'ondes evanescentes 
issues d'un reseau lineaire soit negligcable vue des reseaux voisins. 

En fixant une valeur de 6, la distribution des reseaux lineaires se resume a calculer le nouvel 
espacement le long de I'axe j/a, 

D,^D + {0,5- 5,){D - b), (s = 2, 3, S - 1), 

oil Sj € [0, 1] est une variable aleatoire obeissant une loi de probabilite uniforme. Cette variable 
est propre a chaque reseau lineaire S. A titre d'exemple, la Fig. 7 illustre la repartition des 
reseaux lineaires pour 6 = 0, 9D et 6 = 0, 6D et pour une serie de tirages aleatoires j de 1 a 
10. Pour une bonne lisibilite de la representation seulement 11 reseaux sont illustres. Lorsque 
b = 0, 9D, le deplacement d'un reseau lineaire autour de sa position initiale est inferieur a 0, 5a, 
tandis que le cas 6 = 0, 6D correspond a un deplacement des reseaux lineaires inferieur a 2a 
autour de leur position initiale. 

La Fig. 8 presente les coefficients de reflexion en energie \TZ\^ obtenus pour deux distances 
limites d'approche : 6 = 0,9D et 6 = 0,61?. Ces courbes sont comparees avec celle obtenue pour 
b = D, c'est-a-dire en situation periodique. Nous avons adopte une echelle logarithmique pour 
\TZ\^ afin de mieux juger cette comparaison. Le cas 6 = 0,9a correspond a un faible desordre 
par rapport a la situation periodique. Pour ce faible desordre, la presence de bandes passantes 
et d'arret est toujours observable. Lorsque I'on diminue la distance limite d'approche, 6 = 0,6D, 
la situation physique est differente. Meme si la periodicite est globalement preservee, la Fig. 8 
affiche une diminution des bandes passantes. 

3.2 Dispersion et attenuation des ondes se propageant dans le reseau 
infini 

Le reseau global de la Fig. 5 est maintenant suppose etre constitue d'un nombre infini de 
reseaux lineaires infinis. Ces derniers sont tous identiques et repartis regulierement avec une 
periode D le long de I'axe 2/3. Les fissures de taille 2a ont leur centre localise aux positions 
yi — jd, ?/3 = £D {j,£ — 0,±1,...). II est bien connu de la theorie de Floquet (periodicite ID) 
ou de Bloch (periodicite 2D et 3D) [T^l US] que les equations differentielles avec des coefficients 
periodiques admettent des solutions de la forme : 

U2 (ys) = w {ya) e''^y\ ou w (y^) = w {y^ + D) . (48) 

La grandeur 7 a le caractere d'un nombre d'onde effectif, caracterisant la propagation des ondes 
dans la direction 2/3, dans un milieu elastique homogene equivalent au reseau bidimensionnel 
periodique infini. On dit que 7 est le nombre d'onde de Floquet. Ce nombre d'onde est a valeurs 
complexes et pent s'ecrire : 

7 = 7' + 17", (49) 

oil la partie reelle 7' = lo/c' est la composante suivant 2/3 du vecteur d'onde 7 des ondes se propa- 
geant avec la vitesse c' dans ledit milieu effectif, et la partie imaginaire 7" est leur attenuation. 

II est relativement facile de calculer 7. Achenbach et Li [9j ont presente une methode de calcul 
du nombre d'onde 7 (voir aussi [H]). La derivation formelle de I'equation de dispersion est basee 
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sur la periodicite du solide dans la direction ^3. Une equation de dispersion qui relie 7-D et kxD 
peut s'ecrire sous la forme [5] : 

2T cos (7!?) = (T^ - i?2 - 1) e'''^'" + 2 cos (/ctI?) , (50) 

oil T et i? sont les coefficients de reflexion et de transmission du mode fondamental d'un reseau 
lineaire infini. 

L'existence de solutions reelles et complexes du nombre d'onde de Floquet peut etre in- 
terpretee en termes desdites bandes passantes et d'arret en frequence. Ce sont, respective- 
ment, des parties du spectre 011 les ondes harmoniques sont propagees ou attenuees, et des 
caracteristiques des ondes se propageant dans des milieux multicouches [12]. Sur les bandes pas- 
santes, le nombre d'onde de Floquet est reel, et les ondes sont propagees essentiellement a travers 
le solide non-attenuees. Les bandes d'arret sont des regions oii le spectre est imaginaire pur ou 
complexe. Ici, en vertu de la solution de Floquet (^5)) . les ondes sont attenuees. Physiquement, 
I'attenuation est due a la reflexion progressive d'une onde se dirigeant dans le sens positif aux 
reseaux lineaires successifs du reseau bidimensionnel. 

3.2.1 Etude numerique 

Des courbes dc dispersion associees au reseau infini sont comparees avec les resultats obtcnus 
pour un reseau fini. Cette comparaison n'est valable que si le reseau fini contient suffisamment 
de reseaux lineaires. Avec S — 15, cette condition est satisfaite. Ceci est justifie par le fait que 
les coefficients de transmission (ou de reflexion) revelent les memes bandes interdites que les 
courbes de dispersion (ou d'attenuation) et que celles-ci restent inchangees pour un nombre S 
plus important de reseaux lineaires. 

Pour d = 2, 5 et £) = 4 (la premiere frequence de coupure correspond alors a kxD = Att), la 
Fig. 9 compare les courbes de dispersion (et d'attenuation) du reseau de fissures avec le coefficient 
de transmission en energie du reseau fini. Les courbes de la Fig. 9a font apparaitre une serie des 
bandes d'arret pour lesquelles 7' = ou 7' = n/D. En termes de cristallographie, 7' = corres- 
pond au centre de la premiere zone de Brillouin, et 7' = tt/D a un vecteur d'onde sur la frontiere 
de cette zone. Lorsque 7' — tt/D, un regime d'ondes stationnaires s'etablit. Dans une bande 
d'arret, les ondes se propageant dans le reseau sont attenuees. Inversement, dans les bandes pas- 
santes, les ondes sont dispersives et se propagent dans le reseau sans s'attenuer (7" — 0). Notons 
que les branches imaginaires du spectre d'attenuation doivent tirer leur origine des points sur les 
branches reelles du spectre de dispersion oia la pente est nulle fl8'. La premiere bande d'arret est 
localisee autour d'une frequence telle que la distance D coincide avec une demi-longueur d'onde 
de I'onde incidente : kxD = n. Les bandes d'arret suivantes sont localisees a des multiples entiers 
de cette frequence. La relation ci-avant a ete etablie par Bragg en cristallographie. II s'agit de 
la condition d'interference constructive entre deux ondes reflechies par deux reseaux lineaires 
(ou plans reticulaires) distants de D. Pour un grand nombre de plans reticulaires, une grande 
quantite d'ondes interferent et donnent naissance aux bandes interdites observees. Les largeurs 
des bandes d'arret sont liees a la densite des fissures dans les reseaux lineaires : plus le rapport 
d est petit, plus les bandes d'arret s'elargissent et done plus les bandes passantes sont reduites. 
Pour des reseaux tres denses, les bandes passantes disparaissent completement car chaque reseau 
lineaire reflechit beaucoup plus qu'il ne transmet. Les ondes transmises sont attenuees lors du 
franchissement de chaque reseau lineaire et sont quasi inexistantes apres quelques periodes D 
parcourues dans le reseau bidimensionnel. 

Dans le meme domaine d'etude, la Fig. 9b presente revolution du coefficient de transmission 
en energie |T|^ d'un reseau fini forme de 15 rangees des fissures. A chaque bande d'arret observee 
sur le spectre de dispersion et d'attenuation. Fig. 9a, coincide ici une transmission nulle. Les 
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bandes interdites sont effectivement designees par le coefficient de transmission et des calculs 
menes pour un plus grand nombre de reseaux lineaires n'ont revele aucune modification quant a 
leur forme. Ceci confirme que, pour S* = 15, le reseau global comporte suffisamment de reseaux 
lineaires pour que la propagation d'ondes puisse etre decrite par la theorie de Floquet. Par 
rapport aux courbes de dispersion et d'attenuation, le tracee du coefficient de transmission fournit 
un resultat nouveau lie a la dimension finie du reseau dans la direction : le coefficient de 
transmission presente, au niveau des bandes passantes, des oscillations aperiodiques de nature 
intcrfcrcntielle. Ces oscillations, de type Fabry-Perot, resultent des interferences entre les ondes 
reflechies par les deux interfaces dudit milieu effectif fini. Ces interferences n'ont pas lieu dans les 
bandes d'arret en raison des fortes attenuations qui les caracterisent. Au contraire, les ondes ne 
sont pas attenuees dans les bandes passantes et peuvent done se propager a longue distance dans 
toute I'epaisseur du milieu effectif. Leur forte dispersion explique I'aperiodicite des oscillations 
observees. Si Ton tragait le coefficient de transmission en fonction de 7'-D, les oscillations se 
periodiseraient avec une periode inversement proportionnelle a I'epaisseur h du milieu effectif. 
Les positions des pics d'interference sont en effet donnees par 7' ft, = jtt, avec j — 0, 1, 2, .... 

Conclusion 

Le cadre general dans lequel s'inscrivait ce travail etait la diffusion multiple des ondes 
elastiques antiplanes dans les milieux solides elastiques. L'objectif etait I'etude de la diffusion mul- 
tiple par une multitude de fissures antiplanes periodiquement distribuees. Au cours de la premiere 
partie de ce travail, I'etude de la diffusion par une fissure antiplane a permis de developper une 
equation integrale singuliere de type Cauchy. La forme particuliere de cette equation integrale se 
prete a la mise en oeuvre d'une methode de resolution numcrique adaptee a la determination du 
saut de deplacement a travers la fissure. En employant la formule de representation d'une fissure 
plate et la condition de periodicite d'une rangee de fissures situees sur la meme droite antiplane, 
une equation integrale de frontieres a ete obtenue pour le saut de deplacement a travers une 
fissure dite de reference. Des resultats numeriques pour les coefficients de reflexion et de trans- 
mission out ete presentes en fonction de I'espacement entre fissures, de la frequence d'excitation 
et de Tangle d'incidence. La seconde partie de ce travail avait trait aux reseaux bidimensionnels 
periodiques, constitues d'un empilement d'un nombre fini de reseaux lineaires infinis de fissures 
en milieu elastique. La methode de calcul de la diffusion par un reseau bidimensionnel repose sur 
la decomposition du reseau bidimensionnel en un nombre fini de reseaux lineaires periodiques 
et infinis qui forment alors des plans reticulaires. La diffusion par chacun des reseaux lineaires 
infinis est determinee par un calcul exact de diffusion multiple. La diffusion plan apres plan est 
calculee ensuite par une methode iterative ou a I'aide de la theorie de Floquet. Conformement 
aux theories des milieux periodiques, la presence de bandes interdites et permises dans le spectre 
des frequences a ete observee. Les effets dus a une perturbation de la periodicite ont egalement 
ete analyses. Cette perturbation consistait a introduire un defaut aleatoire sur la position de 
chaque reseau lineaire. L'information principale qui en ressort est que les bandes d'arret sont 
tres sensibles a ce type de perturbation, disparaissant assez tot pour un faible desordre introduit. 

La theorie developpee ici pourra faire I'objet de developpements ulterieurs et etre appliquee 
a un reseau bidimensionnel avec tout type d'aperiodicite, telle qu'une variation des tallies des 
fissures. En autre, en distribuant de fagon aleatoire les reseaux lineaires, le calcul exact des champs 
moyens reflechis et transmis par le milieu aleatoire (les moyennes etant calculees en realisant un 
grand nombre de tirages aleatoires sur les positions des reseaux lineaires) pourra etre compare 
avec des theories dites de milieu effectif (e.g. [19]). Ces theories probabilistes fournissent les 
caracteristiques de I'onde plane, dite coherente, qui se propage dans le milieu homogene equivalent 
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(vu de I'onde coherente) au milieu aleatoire. 
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Fig. 1 : Diffusion d'une onde antiplane par un reseau periodique lineaire de fissures. 
Fig. 1 : Scattering of an antiplane wave by a linear periodic array of cracks. 
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Fig. 2 : Coefficient de refiexion |i?o| versus fa frequence adimensionneiie ui. Effet de fa distance 
d. a) Incidence normaie {9q = 0). b) Incidence obiique {9q = 7r/4). 

Fig. 2 : Reflection coefficient \Rq\ versus the dimensionless frequency uj. Effect of the distance d. 
a) Normal incidence (6q — 0). b) Oblique incidence (Oq — 7r/4j. 
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Fig. 3 : CoefScient de reflexion |i?o| versus la frequence adimensionnelle uj pour une distance 
d — 20 et pour deux valeurs de Tangle d'incidence. Comparaison avec le coefficient de reflexion 
d'une fissure seule. 

Fig. 3 : Reflection coefficient \Ro\ versus the dimensionless frequency uj for a distance d — 20 
and for two values of the incidence angle. Comparison with the reflection coefficient of a single 
crack. 
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Fig. 4 : Coefficient de reflexion |i?o| versus la distance d pour = 0. Influence de la frequence 
adimensionnelle uj. 

Fig. 4 : Reflection coefficient \Ro\ versus the distance d for — 0. Influence of the dimensionless 
frequency a). 
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Fig. 5 : Diffusion d'une onde antiplane par un reseau periodique bidimensionnel de fissures. 
Fig. 5 : Scattering of an antiplane wave by a two-dimensional periodic array of cracks. 
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Fig. 6 : Coefficient de transmission en energie \T\ d'un reseau bidimensionnel de fissures 
excite sous incidence normale avec d = 5 et D = 1. 

Fig. 6 : Energy transmission coefficient \T'f' of a two-dimensional array of cracks excited at 
normal incidence with d = b and D = 1. 
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Fig. 7 : Repartition des reseaux lineaires pour deux distances limites d'approche, 6 = 0, 91? 
et 6 = 0, 6D, et pour 10 tirages aleatoires. 

Fig. 7 : Distribution of linear arrays for two exclusion distances, b = 0.9-D and b = O.QD, and 
for 10 random realisations. 
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Fig. 8 : CoefScients de reflexion en energie \TZ\ d'un reseau bidimensionnel de fissures avec 
d = 5 et D = 1, obtenus pour deux distances limites d'approclie : a) b — 0, 9D et b) 6 = 0, 6D. 
Ces courbes sont comparees avec celle obtenue pour b — D, c'est-a-dire en situation periodique 
(trait fin). 

Fig. 8 : Energy reflection coejficients \Tl\ of a two-dimensional array of cracks with d = 5 and 
D = 1, obtained for two exclusion distances : a) b — and b) b = 0.6D. These curves are 

compared with the one obtained for b — D, i.e. in periodic arrangement (fine line). 
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Fig. 9 : a) Courbes de dispersion (trait fin) et d'attenuation (trait gras) d'un reseau bidimen- 
sionnel infini de fissures avec d = 2, 5 et _D = 4. b) Coefficient de transmission en energie |T| 
d'un reseau bidimensionnel fini de fissures. 

Fig. 9 : a) Dispersion (thin line) and attenuation (solid line) curves of an infinite two-dimensional 
array of cracks with d = 2.5 and D — A. b) Energy transmission coefficient |Tp of a finite two- 
dimensional array of cracks. 
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